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Abstract. The understanding of nonlinear PT-symmetric quantum systems, arising 
for example in the theory of Bose-Einstein condensates in PT-symmetric potentials, 
is widely based on numerical investigations, and little is known about generic features 
induced by the interplay of PT-symmetry and nonlinearity. To gain deeper insights 
it is important to have analytically solvable toy-models at hand. In the present paper 
the stationary states of a simple toy-model of a PT-symmetric system are investigated. 
The model can be interpreted as a simple description of a Bose-Einstein condensate in 
a PT-symmetric double well trap in a two-mode approximation. The eigenvalues and 
eigenstates of the system can be explicitly calculated in a straight forward manner; 
the resulting structures resemble those that have recently been found numerically for a 
more realistic PT-symmetric double delta potential. In addition, a continuation of the 
system is introduced that allows an interpretation in terms of a simple linear matrix 
model. 



PACS numbers: 03.65.Ge, 03.75.Hh 
1. Introduction 

Recent years have seen an increasing interest in non-Hermitian, and in particular PT- 
symmetric quantum theories [TJ [2], which concern quantum systems with complex 
Hamiltonians that nevertheless can have purely real eigenvalues. PT-symmetric 
quantum theories with a purely real spectrum (the case of so-called unbroken PT 
symmetry), have been suggested as generalisations of Hermitian quantum theory on 
a fundamental level [3j. Considerable work has been devoted to the introduction of a 
theoretical framework that leads to unitary time evolutions and physically meaningful 
concepts of observables in these systems (see, e.g., [3, H] and references therein). The 
focus of research in PT-symmetric quantum theories has shifted recently, with the 
proposal and development of experimental applications, where PT-symmetric quantum 
theories arise as effective non-Hermitian descriptions of certain open quantum systems, 
that is, systems in the presence of absorption and gain (see, e.g. [H] and references 
therein). There has been considerable progress in this area, and the first PT- 
symmetric systems have been experimentally realised in optical waveguide structures 
[6j [7] , microwave cavities [8] and electronic circuits [9 J . While there have been a number 
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of proposals for the realisation of PT-symmetric theories in bona fide quantum systems, 
such as Bose-Einstein condensates (BECs) in optical lattices []JJ [TT] [121 E31 EE2 H3 [16] , 
so far these models have not been experimentally realised. 

The consideration of BECs adds an additional interesting feature to the standard 
models of PT-symmetric quantum systems, as the effective description of BECs via a 
mean-field theory is typically nonlinear. Other fields in which nonlinear PT-symmetric 
problems have been investigated are classical wave systems [TT], and complex nonlinear 
Schrodinger dynamics arising in the context of nonlinear optics [TU HE] . The interest in 
the field is still relatively new, and we are only beginning to understand the peculiarities 
of nonlinear versions of PT-symmetric systems. One main hinderance is the lack 
of analytically tractable systems and toy models in which generic properties can be 
analysed and investigated. 

In the present paper we present an analytically solvable nonlinear PT symmetric 
eigenvalue problem that models stationary states of a Bose-Einstein condensate in a PT- 
symmetric two-mode potential. This is the stationary version of a model introduced in 
[TOj, [TTj as the mean- field limit of a PT-symmetric two-mode many-boson system [T3"1, [T9]. 
In [JtJJ [TT] the dynamical properties have been investigated in detail. Here we extend 
these results to give analytic expressions for the corresponding nonlinear eigenvalues, and 
analyse the breakdown of PT-symmetry in this nonlinear model. Eigenvalue structures 
similar to the ones obtained here have recently been found numerically in a model of 
a Bose-Einstein condensate in a double-delta potential [16J. Similar models have also 
been investigated in the context of two coupled nonlinear wave guides [TS] . 

A typical feature of nonlinearity is that it can lead to additional stationary states. 
In the model investigated here, there are up to four eigenstates for arbitrary parameter 
values [TT], while there are only two eigenstates in the linear limit. Allowing for suitable 
complex extensions of the system, four eigenstates can be found for almost all parameter 
values, with the exception of isolated exceptional points at which two stationary states 
coalesce [TJ [20] . These four eigenvalues can be interpreted as those of a linear 4x4 
matrix, which we give explicitly. We furthermore show that there is also a direct relation 
between the stationary states of the nonlinear 2x2 and the linear 4x4 system. 

The paper is organised as follows: In Section [2] we introduce the model system and 
contrast it to similar models considered in the literature. In Section [3] we investigate the 
corresponding time-independent system, provide analytic expressions for the eigenstates 
and eigenvalues, and investigate the stability of the stationary solutions. In Section [4] 
we introduce an analytic continuation of the model that has four eigenstates for all 
parameter values and show that it can be interpreted in terms of a linear matrix model. 
We conclude in Section [5] with a summary and outlook. 

2. The model system 

For low temperatures, the dynamics of a BEC in a real- valued potential V(x) can be 
described in a mean-field approach by the following nonlinear Schrodinger equation 
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(NLSE) of Gross-Pitaevskii type [2T] 

+ v(x) +9 mx,t)\^ t(x,t) = m 9 ^. (i) 

The coefficient of the nonlinear term is given by g = 47rh 2 a/m, where a is the s-wave 
scattering length. 

In the case of a double well potential for low energies it is often appropriate to 
approximate the wave function by a two-mode model by writing 

where £i ; 2 are localised in the left or the right side of the double well, respectively (see, 
e.g., [22J). Neglecting the nonlinear coupling of the two modes, one obtains the two-state 
model 




9\i> 



|2 





H<MM 2 )\ r. 1 . (2) 



with a Hamiltonian matrix H(\ipi\ 2 , \ip2\ 2 ) that depends on the population of the two 
modes. The parameter v is related to the height of the barrier between the two wells, 
which are assumed to be symmetric here. Note that the same Hamiltonian also appears 
in a variety of other physical systems, including nonlinear optics and molecular dynamics 
[231 Ell- 
in the context of BECs, the Gross-Pitaevskii equation is an effective single particle 
approximation of a many-particle quantum system of cold atoms. Although it has 
been introduced formally only for real potentials, recently there has been considerable 
interest in complex generalisations of Gross-Pitaevskii equations [251 EE ETJ EH [161 H2] • 
Considering in particular potentials whose real part is symmetric and whose imaginary 
part anti-symmetric, V(—x) = V*(x), one obtains PT-symmetric versions of the Gross- 
Pitaevskii equationlfj In this spirit, a BEC in a PT-symmetric two-mode potential is 
described by a nonlinear Hamiltonian of the form 



h={ 9 ^ 2 -^ u ; . ]. ( 3) 

1 <# 2 | 2 + i7 1 



This system is PT symmetric with respect to a parity operator P that interchanges the 
two basis states, P : ipi -H- ip2, and the standard time reversal operator T of complex 
conjugation, T : i — > — i. The physical interpretation of the complex energies is a loss 
of particles in the side with negative imaginary part, and an inflow of particles in the 
side with positive imaginary part. The model can be transformed to a purely lossy non- 
Hermitian model via a constant imaginary energy shift [TUl ITT1 [Bj. The Hamiltonian 

| Strictly speaking, the PT-symmetry of the system is only given, if the wave-function respects the 
symmetry and thus renders the nonlinear part in the equation PT-symmetric. 
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Q also appears in the description of PT symmetric wave guide structures with Kerr 
nonlinearity |18j . 

The Gross-Pitaevskii equation, however, is an effective description of a many- 
particle system, valid if the system is condensed; and it is a priori not clear whether 
this so-called mean-field approximation commutes with the complexification of the 
system. The mean-field approximation of non-Hermitian many-particle systems is 
closely related to the classical limit of non-Hermitian single-particle systems, which 
itself is a topic of recent research [29j EQJ El] . In the investigation of example systems 
it has been shown that these two limits do not commute in general. In [TOT ITT] the 
mean-field approximation of a two-mode many-particle system with complex energies 
has been introduced, and it has been found that due to the non- conservation of the 
total probability in the presence of complex energies^] the nonlinear term in the resulting 
Gross-Pitaevskii equation is modified. Instead of the system ^ one finds the dynamical 
equation 



(4) 



Since the norm n = \ipi\ 2 + IV^I 2 is time dependent for non-vanishing 7, the dynamics 
resulting from Q is fundamentally different from the dynamics generated by complex 
Gross-Pitaevskii equations of the form (|3]). Details of the dynamical features of Q can 
be found in pjjj [TT] , whereas the dynamics of ^ is discussed in p3] . Obviously Q and 
(|3j) coincide in both the linear and the Hermitian limit. Further, the time-independent 
versions of these two complex nonlinear Schrodinger equations only differ by a constant 
scaling of the nonlinearity parameter. Thus, when investigating stationary states the 
systems Q and ^ are equivalent. 

3. Stationary states and nonlinear eigenvalues 

In this section we analyse the stationary solutions of the form ip(t) = e~ iflt <j) for 
the nonlinear complex Schrodinger equation Q, and the corresponding generalised 
nonlinear eigenvalues /1. Assuming to be normalised to unity this leads to the time 
independent complex nonlinear Schrodinger equation 



-17 + #101 



| 2 V 

i7 + #|02| 2 





(5) 



with chemical potential /x, which can become complex in the present case. To simplify 
the calculation we introduce the constant energy shift f(|0i| 2 + |02| 2 ), leading to the 
eigenvalue equation 



-i 7 + c(|0i| 2 - |0 2 | 2 ) v 

v i7-c(|0i| 2 - 




>2 



§ Even if the system is in the unbroken PT-symmetric phase the dynamics of an arbitrary initial state 
that is not an eigenstate does not conserve the norm. 



Figure 1. Real (left) and imaginary (right) part of the complex nonlinear eigenvalues 
Q as a function of 7 for v = 1 and c = (top) and c = 0.5 (bottom). 



with c = g/2. 

Note that in the nonlinear case the eigenvalues, that is, the chemical potentials, 
generally differ from the energy expectation values in the stationary states. The 
stationary energy values of the system Q have been investigated in [TTJ. Closely 
related eigenvalue problems have been considered in [261 \27\ [12] . We find that for a 
given stationary solution <ft the corresponding eigenvalue can be expressed as 

fj, = —1JK + ^(0102 + 0200 + CK - 2 - (?) 

The energy in this state, on the other hand, is given by 

E = —i^K + v((f>l4>2 + 020l) + 7^' ( 8 ) 

with k = |0i| 2 — |02| 2 - There are up to four stationary solutions, which can be seen as 
follows. The population imbalance k = \4>i\ 2 — \4>2\ 2 for the possible stationary states is 
given by the real roots of the fourth order polynomial [HI [22] 

(c 2 + 7 V+(^-c 2 -7V = 0, (9) 

which has the four solutions n — 0, 0, ± ^1 — e2 + 7 2 ■ Parameterising the eigenstates as 

0i = /^e--, fc = yiE«er«, (10) 

we find from the eigenvalue equation for k = 

sin(2g) = 1 (11) 
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which leads to two possible values of q. For the eigenvalues k = ±y 1 — e2 ^ 2 , on the 
other hand, there is only one value of q fulfilling 

c 7 
cos (2q) = = and sin (2q) = =. (12) 

V c 2 + 7 2 a/c 2 + 7 2 

In summary, we have the four solutions 

r (0,| arcsin g)) 
(«,?) = ^ (13) 

However, only those pairs k, q which are real valued correspond to actual solutions of the 
time independent Schrodinger equation That is, we have two stationary solutions 
for I7I > \v I and 7 2 < v 2 — c 2 , and four stationary states in the remaining region of the 
parameter space (see also figure 7 in [TT]). 

Clearly, in the linear case there are two eigenstates for each set of parameters. These 
are PT-symmetric if they are symmetric or antisymmetric with respect to the two basis 
states, that is, when k = 0. Thus, PT-symmetry is broken at the critical value of the 
non-Hermitian parameter \^ C ru\ = M- The two corresponding eigenvalues, depicted in 
the upper panel in Figure [lj are simply given by fi± = ±\Jv 2 — 7 2 . They are real in the 
region of unbroken PT symmetry, and become complex when the eigenstates break the 
PT symmetry. 

The lower panel of Figure [T] shows the nonlinear eigenvalues for a relatively small 
nonlinearity. We observe a similar behaviour as in the linear case, with two real 
eigenvalues that coalesce in an exceptional point at the critical value I7I = \v\. These 
eigenvalues are independent of the nonlinearity, as they correspond to k = and 
thus the nonlinear term vanishes in the eigenvalue equation (|6j). However, instead of 
turning to a complex conjugate pair after the exceptional point, the eigenvalues and 
the corresponding eigenstates vanish in the nonlinear case. On the other hand, there is 
an additional pair of complex conjugate eigenvalues, which emerge already for smaller 
values of |7|. 

In general, in the nonlinear case c 7^ there can be up to four coexisting eigenstates, 
only two of which respect the PT-symmetry of the system. Figure [2] shows examples 
of the eigenvalues in dependence on 7 for different values of c. Since the imaginary 
part of the eigenvalue is proportional to the population imbalance k, PT symmetric 
solutions correspond to real eigenvalues, as in the linear case. The solutions with k^O 
appear for values of the non-Hermiticity parameter for which 7 2 > v 2 — c 2 . Thus, in the 
presence of nonlinearity the PT-symmetry of the system is broken for smaller values 
of the non-Hermiticity, and above the critical value \c cr u\ = \v\ there is no region of 
unbroken PT symmetry at all. 

In the Hermitian case additional states induced by the nonlinearity are often 
referred to as self-trapping states, because they favour one of the modes, thus breaking 
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Figure 2. Real (top) and imaginary (bottom) part of the complex nonlinear 
eigenvalues (J7]) as a function of 7 for v = 1 and c = 0.9 (left), c = 1 (middle), 
and c = 1.5 (right). 



the symmetry of the linear part of the system. Therefore, one can argue that the non- 
Hermiticity enhances the self trapping effect, by reducing the critical nonlinearity for 
which additional states emerge, however, on cost of the PT-symmetry of the system. 

It is well known that the self-trapping transition in the two-mode Gross-Pitaevskii 
equation with real parameters is accompanied by a restructuring of the dynamical 
stability of the eigenstates: While the self-trapping modes that break the symmetry of 
the system are dynamically stable, the previously stable antisymmetric state becomes 
dynamically unstable after the self-trapping transition. The stability can be determined 
via the Bogoliubov-de-Gennes equation for the stationary states, which arise from the 
linearised effect of a small time-dependent perturbation of the stationary solution of the 
form if)(t) = e~ tfJjt (j) + <50_e~ lwi + 5(f* + e luJ *V However, in the case of complex chemical 
potential, care has to be taken to consider perturbations that are orthogonal to the 
stationary solution (f> j32l [28] . In the present case the resulting Bogoliubov-de-Gennes 
equations reduce to the eigenvalue problem 
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with the 4x4 matrix 
( 



B- 




Bo«-p+cQy_ m ^ JQ cQ 




2 



V 

Here Q denotes the projection operator orthogonal to the stationary solution </>, 

v 

and H Qn i is the complex nonlinear Hamiltonian evaluated at the stationary solution 

n ° nl ~ „. „f\A 12 IJ. I2N 




i7-c(|0i| 2 - 



>1\ 



The stability of a stationary solution is determined by the imaginary part of the 
eigenvalues u: While for lm(u) < the perturbations are decaying or purely oscillating, 
they grow exponentially for lm(u) > 0, in which case the stationary solution is 
dynamically unstable. 

Figure [3] shows the numerically obtained eigenvalues u in the complex plane for 
all stationary states for c = 0.5, v = 1 and three different values of 7. The figure 
in the left panel corresponds to 7 = 0.2, where there are only two eigenstates with 
real chemical potentials, which are both stable. For 7 = 0.7, already in the region 
of broken PT-symmetry, there are two additional eigenstates with complex conjugate 
chemical potentials. The one with a positive imaginary part of the chemical potential 
is dynamically stable, while the other one is unstable. This is related to the fact 
that these states appear as local sink and source of the dynamics [TD1 [EE] . As in the 
Hermitian case, one of the remaining stationary states is now unstable, it is a saddle 
point of the dynamics [TUl E]. It can be easily verified analytically that this state 
changes its dynamical behaviour at the point in parameter space at which the additional 
eigenstates emerge: The eigenvalue equation for the stability matrix B for the states 
(n,q) = (0, I arcsin C 1 )), that correspond to the chemical potentials \i = ±^/v 2 — 7 2 
reduces to 

co A + 4( 7 2 - v 2 ± c^v 2 - 7 2 V 2 = 0. (15) 



Thus, for positive c the upper state becomes dynamically unstable for 7 > \/v 2 — c 2 , 
while the other state stays stable. The situation is reversed for negative values of c, 
in this latter case it is also the lower chemical potential from which the additional 
eigenvalues branch off. This behaviour is in complete analogy to the self-trapping 
transition for 7 = 0. For even larger values of 7 > 1 only the two eigenstates with 
complex chemical potential remain, which do not change their stability, that is one of 
them is stable, the other is unstable. 
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Figure 3. Eigenvalues cu of the stability matrix B corresponding to the different 
stationary states in the complex plane for c = 0.9, v = 1 and 7 = 0.2 (left), 7 = 0.7 
(middle), and 7 = 1.2 (right). 



We have seen that the complex nonlinear eigenvalue problem (|6]) leads to a varying 
number of solutions, depending on the precise parameter values, where eigenvalues may 
appear, vanish and change their stability when parameters are varied. While such 
bifurcations are typical in nonlinear problems, it might lead to useful insights to consider 
an extension of ^ that has four eigenvalues for all parameter values, which coincide 
with the eigenvalues of the original system wherever the latter are defined. We will 
introduce such an extension in the following section. 



4. Analytic extension of the Gross-Pitaevskii equation and an isospectral 
linear system 

The vanishing of solutions in the Gross-Pitaevskii equation is related to the fact that 
the nonlinear term g\ip\ 2 is not analytic [331 US], and the "missing" eigenstates and 
eigenvalues can be recovered by considering analytic extensions of the problem. The idea 
presented in [33] is as follows: Every wave function ip(x) can formally be decomposed in 
terms of two real valued functions A(x) and (p(x) corresponding to the amplitude and 
the phase of ip(x) = A(x)e ltp ^ x \ The non-analytic complex Gross-Pitaevskii equation 
can thus be decomposed into two coupled real differential equations. To extend the 
parameters into the complex domain, rather than considering the original equation, one 
formally considers the equations for A(x) and (p(x), which now lead to complex solutions. 
This is equivalent to replacing |t/>(x)| 2 in the Gross-Pitaevskii equation with the complex 
valued function A 2 (x). In the present case, this leads to the complex Gross-Pitaevskii 
equation 

-i7 + ck v 



V 17 — CK 




H - (16) 



where the eigenstates 02 ) are parameterised according to (10), with arbitrary 



k,5 6 C. When solving the eigenvalue equation we have in fact already encountered the 



additional solutions, which led to complex values of k or q in (13). The corresponding 
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Figure 4. Real (top) and imaginary (bottom) part of the complex continued mean- 
field eigenvalues (17) as a function of 7 for v = 1 and c = 0.5 (left), c = 1 (center) and 
c = 1.5 (right). 



eigenvalues are then given by 

H = (— 'vy + ck)k + vyl — k 2 cos(2g), (17) 



for all solutions (k, q) in (13). Substituting the values of k and q into (17) we find the 
eigenvalues 




V={ _ , _. / ^ 1 (18) 

c 2 +7 2 J- , 

that is, the eigenvalues of the linear problem and two additional nonlinear eigenvalues. 
Thus, the eigenvalues are naturally divided into two separate sets. The resulting 
spectrum, as a function of 7, is depicted in Figure |4j In the limit c — » the two 
nonlinear eigenvalues approach the linear ones. It should be noted that in this limit also 
the nonlinear eigenstates coincide with the linear ones, and thus we recover two copies 
of the linear system. 

It is an interesting question whether spectral features of nonlinear systems can be 
mimicked by suitable linear systems, and which symmetries these linear systems exhibit. 
In the present case the spectrum of the analytically continued nonlinear eigenvalue 



problem (16) can be immediately interpreted as that of an equivalent four-dimensional 
linear system. Since the eigenvalues decouple into two sets, it is not difficult to construct 
a matrix representation of this linear system. It is straightforward to verify that the 
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nonlinear problem (16) is isospectral to the linear 4x4 matrix: 



M 



( -Y) V \ 

v i 7 

-vy + c 

V o ^ i 7 + c/ 



1JV 

c+i-y 



(19) 



This matrix is itself PT-symmetric. However, there is an interesting subtlety to this: 
While the upper block has the same symmetry as the linear two-level system, the lower 
block has additional complex elements on the off-diagonal, and is not PT symmetric 
with respect to the standard T operator defined as complex conjugation alone. Instead, 
the second block is PT symmetric with respect to the P operator exchanging the two 
relevant basis elements, and the T operator defined as Hermitian conjugation. 

Let us finally turn to the relation between the eigenstates of the two systems. On 
a formal level, they are in fact directly related: The eigenstates of the nonlinear two- 



level system ( 16 ) appear in the upper and lower components of the eigenvectors of the 
matrix (19), while the remaining components automatically vanish due to the block 
diagonal form of the system. This, however means, that there is no coupling or overlap 
between the two sets of eigenstates, and in particular the degeneracy between one of 
the eigenvalues of the upper block with the degenerate pair of eigenvalues of the lower 
block at 7 = \Jv 2 — c 2 has to be interpreted as ocuring incidentally in the linear system. 
Since the dynamics of the two systems are obviously different, it is not surprising that 
also the dynamical stability properties of the eigenstates differ. In the linear model it 
can be checked in a straightforward manner that all eigenstates are stable in the region 
of unbroken PT-symmetry. This includes the Hermitian limit 7 — > 0, which means that 
the reorganisation of the stability properties due to the self-trapping transition is not 
reflected in the linear system. On the other hand, once the PT-symmetry is broken, 
only the eigenstate corresponding to the eigenvalue with the largest imaginary part is 
dynamically stable in the linear system, leaving the remaining three states unstable. 
This is connected to the fact that the eigenstate corresponding to the eigenvalue with 
the largest imaginary part acts as a global sink of the dynamics. Thus, the equivalence 
between the linear and the original nonlinear complex eigenvalue problem clearly has 
limitations, in particular in regard to the related dynamical properties. It is, however, 
an interesting question whether there might be a hidden deeper relation between the 
models, for example in a sense of constraint quantum mechanics |34j . 

The final interpretation of the isospectral linear operator for the nonlinear 
eigenvalue problem is an open problem, the relevance of which will only become apparent 
in the study of more general systems. A deeper understanding of the relation between 
the two systems might lead to valuable insights into the nature of both real and complex 
nonlinear eigenvalue problems, and the nature of complex symmetries. The investigation 
of these issues, however, goes beyond the scope of the present paper. 
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5. Summary and Outlook 

We have presented an analytically solvable paradigm system for complexified nonlinear 
quantum systems. The model can be interpreted as a simple description of a BEC in a 
PT symmetric double well potential. In particular, we have identified and analysed the 
eigenvalues and corresponding eigenstates of the system and their stability properties. 
Similar to the Hermitian case, the nonlinearity can lead to additional eigenstates, which 
in the complexified case leads to a breaking of PT-symmetry. This analytically solvable 
model is expected to be valuable in the understanding of more complicated systems, 
which can only be investigated numerically. An example is the model of a BEC in a 
non-Hermitian double delta potential, numerically investigated in [16J, where similar 
features to those discussed here have been found. 

We have further introduced a continuation of the problem that is isospectral to 
a simple matrix model, whose eigenstates are also related to those of the nonlinear 
system. In the present case this system was directly constructed from the known 
solutions of the nonlinear system. If a more general connection between these isospectral 
models could be uncovered this would shed new light onto the interpretation of nonlinear 
eigenvalue problems. To gain further insights into this connection, it would be useful to 
investigate more general systems, for example, discrete nonlinear Schrodinger equations 
with additional modes or systems with different types of nonlinearities. 

I thank Giinter Wunner for stimulating discussions, and Hugh Jones, Dorje Brody, 
and Hans Jiirgen Korsch for useful comments. I would like to express my gratitude 
to Nimrod Moiseyev and the Chemistry Department of the Technion, Haifa for kind 
hospitality while this manuscript was prepared. I am also grateful for support from the 
Imperial College Junior Research Fellowship scheme. 
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